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Phase transitions waves in atomic chains with double-well potential play a fundamental role in 
materials science, but very little is known about their mathematical properties. In particular, the 
only available results about waves with large amplitudes concern chains with piecewise-quadratic 
pair potential. In this paper we consider perturbations of a bi-quadratic potential and prove that 
the corresponding three-parameter family of waves persists as long as the perturbation is small 
and localised with respect to the strain variable. As a standard Lyapunov-Schmidt reduction 
cannot be used due to the presence of an essential spectrum, we characterise the perturbation 
of the wave as a fixed point of a nonlinear and nonlocal operator and show that this operator is 
contractive in a small ball in a suitable function space. Moreover, we derive a uniqueness result 
for phase transition waves with certain properties and discuss the kinetic relation. 
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1 Introduction 

Many standard models in one-dimensional discrete elasticity describe the motion in atomic chains 
with nearest neighbour interactions. The corresponding equation of motion reads 



where ^ is the interaction potential and Uj denotes the displacement of particle j at time t. 

Of particular importance is the case of non-convex because then (I) provides a simple dy- 
namical model for martensitic phase transitions. In this context, a propagating interface can be 
described by a phase transition wave, which is a travelling wave that moves with subsonic speed and 
is heteroclinic as it connects periodic oscillations in different wells of The interest in such waves is 
also motivated by the quest to derive selection criteria for the naive continuum limit of (1), which is 
the PDE dttu = dx^'{dx)- For non-convex this equation is ill-posed due to its elliptic-hyperbolic 
nature, and one proposal is to select solutions by so-called kinetic relations [AK91, Tru87] derived 
from travelling waves in atomistic models. 

Combining the travelling wave ansatz Uj{t) = U{j — at) with (1) yields the delay-advance- 
differential equation 
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where R{x) := U{x + 1/2) — U{x — 1/2) the (symmetrised) discrete strain profile and AiF{x) : = 
F{x + 1) — 2F{x) + F{x — 1). Periodic and homoclinic traveUing waves have been studied intensively, 
see [FW94, SW97, FP99, Pan05, EPOS, HerlO] and the references therein, but very little is known 
about heterochnic waves. The authors are only aware of [HRIO, Herll], which prove the existence of 
supersonic heterochnic waves, and the small amplitude results from [looOO]. In particular, there seems 
to be no result that provides phase transitions waves with large amplitudes for generic double-well 
potentials. 

Phase transition waves with large amplitudes are only well understood for piecewise quadratic 
potentials, and there exists a rich body of literature on bi-quadratic potentials, starting with [BCSOla, 
BCSOlb, TV05]. For the special case 

$(r) = — \r\ , $'(r) = r — sgn(r) (3) 

the existence of phase transition waves has been established by two of the authors using rigorous 
Fourier methods. In [SZ09] they consider subsonic speeds c sufficiently close to 1, which is the speed 
of sound, and show that (2) admits for each c a two-parameter family of waves. These waves have 
exactly one interface and connect different periodic tail oscillations, see Figure 2 for an illustration. 

In this paper we allow for small perturbations of the potential (3) and show that the three- 
parameter family of phase transition waves from [SZ12] persist provided that the perturbation is 
sufficiently localised with respect to the strain variable r. 

A related problem has been studied in [VailO]. There, a piecewise quadratic family of potentials 
is considered such that the stress-strain relationship is continuous and trilinear, with a small spinodal 
region. Travelling wave solutions arc shown to obey a relation of residuals in the Fourier represen- 
tation, which is then approximately solved numerically. The regularity of the perturbed potential is 
lower than that of the class of perturbations considered here, so strictly speaking the results do not 
overlap. However, in spirit the settings are close and indeed the numerical evidence [VailO, Fig. 4, 
bottom right panel] is in good in agreement with our findings: there is an one-sided asymptotically 
constant solution, and the tail behind the interface oscillates with slightly different amplitude from 
that related to (3). The range of velocities considered in [VailO] is larger than the one studied here. 

Our approach is in essence perturbative and reformulates the travelling wave equation with per- 
turbed potential in terms of a corrector proffie S, i.e., we write R = Rq + S, where Rq is a given 
wave corresponding to the unperturbed potential. The resulting equation for the corrector S can be 
written as 

MS = A^g{S) + r], (4) 

where r/ is a constant of integration and A, A4, G are operators to be identified below. More precisely, 
is a linear integral operator which depends on c and Q a nonlinear superposition operator involving 
Rq. The analysis of (4) is rather delicate since the Fourier symbol of M. has real roots, which implies 
that is an inner point of the continuous spectrum of M. In particular M. is not a Predholm 
operator in the function spaces we considered here, so a standard bifurcation analysis from 6 = via 
a Lyapunov-Schmidt reduction is impossible. 

In our existence proof, we first eliminate the corresponding singularities and derive an appropriate 
solution formula for the linear subproblem. Afterwards we introduce a class of admissible functions 
S and show that the properties of A and Q guarantee that A'^g{S) is compactly supported and 
sufficiently small. These fine properties are illustrated in Fig. 5 and allow us to define a nonlocal and 
nonlinear operator T such that 

MT{S) = A'^GiS) + r]{S) 

holds for all admissible S with some r]{S) € M. This operator T is contractive in some ball of an 
appropriately defined function space, so the existence of phase transitions waves is granted by the 
contraction mapping principle, see Lemma 14. Moreover, the properties of Ai and Q imply that our 
fixed point method for S yields all phase transition waves R that comply with certain requirements, 
see Proposition 17. 
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Our existence result yields - for each c from an interval of subsonic velocities - a genuine two- 
parameter family of solutions to (2) but it is not clear whether all these phase transition waves are 
physically reasonable. In the literature, one often employs selection criteria to single out a unique 
phase transition wave for each speed c. One selection criterion is the causality principle, which in 
our case selects waves with non-oscillatory tails in front of the interface; see [SlcOl, Slc02, TV05] 
and Remark (v) following the Main Theorem 3. These waves can also be observed in numerical 
simulations of atomistic Riemann problems with non-oscillatory initial data [HSZ12]. 

Below we tailor our perturbation method carefully in order to show persistence of the tail oscil- 
lations in front of the interface. In particular, for each small 6 and any given c we obtain exactly one 
wave that propagates towards an asymptotically constant state. The other solutions are oscillatory 
for both X — >■ — oo and x — >■ -|-oo, and satisfy the entropy principle - which is less restrictive than the 
causality principle - as long as the oscillations in front of the interface have smaller amplitude than 
those behind; see [HSZ12] for more details and a discussion of the different versions of Sommerfeld's 
radiation condition. It is not known whether waves with tail oscillations on both sides of the inter- 
face are dynamically stable or can be created by Riemann initial data. As usual, however, one might 
expect that travelling waves resemble dynamical solutions of more complex situations in a temporal 
and spatial window. Candidates arc cascades of moving phase interfaces in chains with multi-well 
potential or other types of macroscopically self-similar waves. 

We also emphasise that phase transition waves satisfy Rankine-Hugoniot conditions for the macro- 
scopic averages of mass, momentum, and total energy [HSZ12], which imply nontrivial restrictions 
between the wave speed and the tail oscillations on both sides of the interface. Although these con- 
ditions do not appear explicitly in our existence proof, they can (at least in principle) be computed 
because the tail oscillations are given by harmonic waves, see again Fig. 2. For general double-well 
potentials, however, it is much harder to evaluate the Rankine-Hugoniot conditions and thus it re- 
mains unclear which tail oscillations can be connected by phase transition waves. Closely related to 
the jump condition for the total energy is the kinetic relation, which specifies the transfer between 
oscillatory and non-oscillatory energy at the interface and determines the configurational force that 
drives the wave. In the final section we discuss how the kinetic relation changes to leading order 
under small perturbations of the potential (3) . 

We now present our main result in greater detail. 

1.1 Overview and main result 

We study an atomic chain with interaction potential 

^5{r) = y -^5{r), *5(0) = 0, 

where ^'^ is a perturbation of ^'q = sgn in a small neighbourhood of 0. The travelling wave equation 
therefore reads 

c^B!' = 1^^{R-%{R)) (5) 

and depends on the parameters c and 5. In order to show that (5) admits solutions for small 5 we 
rely on the following assumptions on see Figure 1 for an illustration. 

Assumption 1. Let (^s)^^^ he a one-parameter family of -potentials such that 

(i) ^'^ coincides with outside the interval {—6, 5), 

(ii) there is a constant C$ independent of 5 such that 

for all r G M. 
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The quantity 



Is :-- 



[^'s{r)-%ir))dr 



plays in important role in our perturbation result as it determines the leading order correction. 
Notice that our assumptions imply 

Is = l J ^ ^'s{r) dr = -|(^>5(+1) - and hence \Is\ < C^S . 

As already mentioned, the case 6 = has been solved in [SZ09]. The main result can be summarised 
as follows. 

Proposition 2 ([SZ09], Proof of Theorem 3.11 and [SZ12], Theorem 1). There exist < cq < 1 
such that for every c G [cq, 1), there exists a two-parameter family of solutions Rq € W^'°^(R) to the 
travelling wave equation (5) with 6 = 0. This family is normalised by Ro{0) = and can be described 
as follows: 

(i) There exists a unique travelling wave Rq such that 

Roix) ft 
Ro{x) — a" ( cos (kcx) — l) — sin {kcx) °° > f~ 

for some constants r^ , kc, a~ , and j3~ depending on c. 

(a) There exists an open neighbourhood Uc of in such that for any (a, /3) G Uc the function 
Rq = Rq + a( cos (kc-) — l) + ^ sin {kc-) is a travelling wave with 

(a) \\Ro\\oo < Do{l - c^y\ 

(b) Ro{x) > ro for x > xq and Ro{x) < — ro for x < —xq , 

(c) R'q{x) > do for \x\ < Xq 

for some constants xq, ro, do, and Dq depending on co- 
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Figure 1: Sketch of and $5 for 5 = (grey) and 6 > (black). Since $0 is symmetric, — 15 is just 
half the energy difference between the two wells of $5. 



The main result of this article can be described as follows. 

Theorem 3. For all ci G (co, 1) there exists 60 > such that for any < 6 < 60, any speed 
Co < c < c\, and any given wave Ro as in Proposition 2 there exists a solution R to (5) with 

R = Ro — Is S, 

where Ig = 0{6) and the corrector S G W^'°°(M) 
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FigTirc 2: Sketch of the waves for 5 = (grey) and S > (black) as provided by our perturbation 
result; the shaded region indicates the spinodal interval [—6, +6], where differs from \I/q. Both 
waves differ by the constant Is = 0{S) and a small corrector S of order 0[6^) , which is oscillatory for 
X < but asymptotically constant as x — >■ +oo. The tail oscillations of both waves do not penetrate 
the spinodal region and are generated by harmonic waves with wave number kc- For each admissible 
6 and c there exists exactly on wave that satisfies the causality principle as it is non-oscillatory for 
X — >■ +00. 

(i) vanishes at x = 0, 

(ii) is non-oscillatory as x +00, i.e., the limit limj;_^+oo S{x) is well defined, 

(Hi) admits harmonic tail oscillations for x — >■ —00, that means there exists constants a_ and d_ 
such that limj;_^_oo S{x) — a-Ro{x + d-) is well defined, 

(iv) is small in the sense of 

ll^lloo = 0(<52), \\S'\\oo = Oi6), \\S"\\oo = 0{l). 

Moreover, the solution R with these properties is unique provided that S is sufficiently small. 

More detailed information about the existence and uniqueness part of our result are given in 
Proposition 15 and Proposition 17, respectively. We further mention: 

(i) Since the travelling wave equation is invariant under 

c — c, R{x) R{—x), 
there exists an analogous result for — 1 < c <C 0. 

(ii) Different choices of c and Rq provide different waves R, see Section 4. 

(iii) The travelling wave equation (5) is, of course, invariant under shifts in x but fixing Rq and S 
at removes neutral directions in the contraction proof. 

(iv) All constants derived below depend on ci and cq but for notational simplicity we do not write 
this dependence explicitly. It remains open whether Sq can be chosen independently of ci. 

(v) The causality principle selects those solutions with Cgr < Cph and Cgr > Cph for all oscillatory 
harmonic modes ahead and behind the interface, respectively, where, Cgr and Cph are the group 
and the phase velocity. For nearest neighbour chains with interaction potential $0 ^'id wave 
speed c sufficiently close to 1, Proposition 2 yields 

Cph = c = a{kc) = k~^U{kc) > Cgr = ^'{kc) 

on both sides of the interface, where i}{k) = 2sin(A;/2) is the dispersion relation [SCC05, 
HSZ12]. The causality principle therefore selects the solution ^0 ^ it is the only wave having 
no tail oscillations ahead of the interface. Since our perturbative approach changes neither the 
wave speed c nor the wave number kc in the oscillatory modes (but only the amplitude behind 
the interface and, of course, the behaviour near the interface), we conclude that Theorem 3 
provides for each 6 and c exactly one wave that complies with the causality principle. 
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(vi) The surprisingly simple leading order effect, that is the addition of —Is to i^o, implies that 
the kinetic relation does not change to order 0{S). Notice, however, that the kinetic relation 
depends on the choice of Ro, cf. [SZ12]. 



This paper is organised as as follows. In Section 2 we reformulate (5) in terms of integral operators 
A and M and show that it is sufficient to prove the existence of waves for the special case Ig = 0. 
Section 3 concerns the existence of correctors S. We first establish an inversion formula for Ai 
which in turn enables us to define an appropriate solution operator C to the affine subproblem 
MS = A?G + T] with given G. Afterwards we investigate the properties of the nonlinear operator Q 
and prove the contractivity of the fixed point operator T. In Section 4 we establish our uniqueness 
result and conclude with a discussion of the kinetic relation in Section 5. 



2 Preliminaries and reformulation of the problem 

In this section we reformulate the travelling wave equation (5) in terms of integral operators and 
show that elementary transformations allow us to assume that I5 = Q holds for all 5 > Q. 

2.1 Reformulation as integral equation 

For our analysis it is convenient to reformulate the problem in terms of the convolution operator A 
and the operator M. defined by 

x+1/2 

{AF) {x):= J F{s)ds, MF := A^F - c^F . 

x-1/2 

The travelling wave equation can then be stated as 

MR = A^^'s{R) + ^Ji. (6) 

Lemma 4. A function R G W^'°°(M) solves the travelling wave equation (5) if and only if there 
exists a constant /x G M such that {R, n) solves (6) . 

Proof. By definition of „4, we have -^A? = Ai. Equation (5) is therefore, and due to the definition 
of M, equivalent to 

{MR)" = P", P := A'^^'siR) . (7) 

The implication (6) (5) now follows immediately. Towards the reversed statement, we integrate 
(7)1 twice with respect to x and obtain MR = P + Xx + fi, where A and n denote constants of 
integration. The condition R £ L°°(M) implies MR, %{R), P G L°°(M), and we conclude that 
A = 0. □ 

2.2 Properties of the operators A and M. 

Some of our arguments rely on Fourier transform, which we normalise as follows 

F{k) = ^ [ e^'=^F(x)dx, F{x) = ^ f e-^^^F(fc)dA:. 
V27r Jm v27r 7m 

Using standard techniques for the Fourier transform in the space of tempered distributions we readily 
verify the following assertions. 
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Figure 3: The the real roots of the symbol function m are the solutions to \a{k)\ = \c\. 



Remark 5. The operators A and M. diagonalise in Fourier space and have symbols 



a{k) = 



sin(/c/2) 
k/2 



and m{k) = a{k) — , 



respectively. In particular, we have 

Mcos{kc-) = 0, Msm{kc-) = 0, Ml = 1 - , 
for any real root kc of m, and 

F G span { cos {kc-), sin [kc-) : m{kc) = 0, kc> O} 
for any tempered distribution F with MF = 0. 

The set of real roots of m depends strongly on the value of c, see Figure 3. In what follows we 
only deal with positive and near sonic speed c, that means c ^ 1, for which m has two simple real 
roots. 

We next summarise further properties of the operator A and recall that the Sobolev space W^"^(M) 
is for any 1 < p < oo continuously embedded into BC(M). 

Lemma 6. For any 1 < p < og we have A: Lp(M) W1'P(M) c BC(R) with 



\AF\\p < \\F\\p , \\{AFy\\p < 2\\F\\p , \\AF\\^ < \\F\\p 



(8) 



for all F G IP(R), where {AF)' = VF := F[- + 5) - F[- - Moreover, suppF C [xi, X2] implies 

SUpp AF C [xi - |, X2 + 5]. 



Proof. Let 1 < p < 00 and F G L*'(M). The definition of A ensures that AF has in fact the weak 
derivative V-F, and this implies the estimate (8)2 via ||VF||p < 2||F||p. Using Holder's inequality we 
find 

|(^F)(x)|^< r'^'\F{s)fds 
Jx-1/2 

and integration with respect to x yields (8)1. We also infer that |(.AF)(a;)| < \\F\\p holds for all 
X G M, and this gives (8)3. Finally, the arguments for p = 00 are similar and the claimed relation 
between supp F and supp AF is a direct consequence of the definition of A. □ 

2.3 Transformation to the special case Is = 

The key observation that traces the general case Ig ^ back to the special case = is that any 
shift in ^'g can be compensated by adding a constant to R. 
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Lemma 7. The family (^^)j>o defined by 

satisfies Assumption 1 with constant C-^ = C^(l + C^) as well as 

i~g = l [ ^'-Ar) - *o(r) dr = for all S > 0. 

Moreover, each solution {R, jl) to the modified travelling wave equation 

MR = A'^^iR) + A (9) 

defines a solution {R, /j,) to (6) via R = R — Is and fj, = jl — (^c^ — 1) Is a,nd vice versa. 

Proof. Due to \Is\ < C-q,S and our definitions we find ^~(r) = ^'o(r) at least for all r with \r\ > S, as 
well as 



5^ 

We also have 



*Ur) <C^<C^, mr) 



C-q, (7,1, 1 + C^, C^j 

< = —r- z = — ^ for all r G 

~ 5 5 l + C^ s 



h = \l - h) - %{r)) dr=U {%ir) - %{r + Is)) dr 

= I [ (*5 W - *o W) dr + i / (*'o(r) - %{r + Is)) dr = Is - Is = . 

JR JR 

Finally, the equivalence of (6) and (9) is obvious. □ 

3 Existence of phase transition waves 

In this section, we show that each phase transition wave for persists under the perturbation 
^'o ^'j, provided that S is sufficiently small. To this end we proceed as follows. 

(i) We fix c G [co, ci] with < cq < ci < 1 as in Proposition 2 and Theorem 3. Then there exists a 
unique solution fcc > to a{kc) = c, and this implies m{±kc) = 0, m'{±kc) 7^ 0, and m{k) 7^ 
for A; 7^ zb/cc- All constants derived below can be chosen independently of c but are allowed to 
depend on cq and ci. 

(ii) Thanks to Proposition 2 and Lemma 4, we fix (Rq, fio) from the two-parameter family of 
solutions to the integrated travelling wave equation (6) loi 6 = and given c. Recall that Rq 
is normalised by Ro(0) = 0. 

(iii) In view of Lemma 7, we assume that Is = holds for all 6 > 0. To avoid unnecessary 
technicalities we also assume from now on that 5 is sufficiently small. 

In order to find a solution (R, ji) to the integrated travelling wave equation (6) for (5 > 0, we further 
make the ansatz 

R = Rq + S , 11 = 110 + r], 
and seek correctors {S, rj) such that 

MS = A^G + ri, G = g{S). (10) 
Here, the nonlinear operator Q is defined by 

g{S){x) = ^'s{Ro{x) + S{x)) - ^'o{Ro{x)) . (11) 
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In order to identify a natural ansatz space X for S, we first remark that the smoothing properties of 
A, see Lemma 6, imply S G W^'°^(M). Notice, however, that R = Rq+S is in general more regular due 
to the smoothness of ^s- More precisely, (6) combined with '^g G C'^(M) yields R G 0*^+^(1^). We also 
impose the normalisation condition S{0) = in order to eliminate the non-uniqueness that results 
from the shift invariance of the travelling wave equation (6). In fact, without this constraint any 
corrector S provides a whole family of other possible correctors via S = S{- + xq) + Ro{- + xq) — Rq 
with Xq = 0{6^). 

A key property of our existence and uniqueness result is that R has only harmonic tail oscillations 

with wave number kc and that both R and Rq share the same tail oscillations for x — )■ +oo. The 
corrector S is therefore non-oscillatory in the sense that S{x) converges as a:; — >■ -t-oo to some well- 
defined limit a. In summary, we seek solutions (S, r/) to (10) with 5 G X and r/ G M, where 

X := 1^ G W2'°°(M) : S{0) = , a = ^Hrn^ S{x) exists } 
is a closed subspace of W^'°°(M) and hence a Banach space. 
3.1 Inversion formula for A4 

Our first task is to construct for given G a solution {S, rf) to the afHne equation (10)i. In a preparatory 
step, we therefore study the solvability of the equation 

MF = Q (12) 

using the Fourier transform for tempered distributions, where Q G L°^(R) is some given function. This 
problem is not trivial because the symbol function m has two simple roots at zt/cc, or, equivalently, 
because is an element of the continuous spectrum of M. corresponding to a two-dimensional space 
of generalised eigenfunctions. We are therefore confronted with the following two issues in Fourier 
space: 

(i) F is uniquely determined only up to elements from the space 

span{ (5_fc^(/c) , 5j^.kS^) }> 
which contains the Fourier transforms of all bounded kernel functions of M. 

(ii) F exhibits - at least for generic Q with (^(ibfcc) 7^ - two poles at ifcc and is hence not 
Lebesgue integrable in the vicinity of ifcc. In particular, the dual pairing between F and a 
Schwartz function is defined in the sense of Cauchy principal values only. 

The non-uniqueness is actually an advantage because it allows us to select solutions with particular 
properties; see the proof of Lemma 10, where we add an appropriately chosen kernel function to 
ensure non-oscillatory behaviour for x — > -|-oo. Concerning the non-integrable poles at ibfcc, we split 
F into a two-dimensional singular part and a remaining regular part, and show that any solution F 
to (12) belongs to some Lebesgue space provided that Q is sufficiently regular. 



graph of Yi graph oi MYi graph of Yi 




Figure 4: Properties of Yi (grey) and Y2 (black). 
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As illustrated in Figure 4, we introduce two functions Yi,Y2 G L°° (M) with 



v27r V27r 
^i(^) •= //, N COS (A;ca;)sgn(x) , Y2{x) := sin (A;cx)sgn(x) , 

m'[kc) m'[Kc) 

and verify by direct computations the following assertions. 

Remark 8. We have 

(i) MYi G L°°(M) with supp MYi C [-1, 1] , 

("m; mfj e L2(M) n BC^(M) mi/i lim m(fc)Yi(A;) = ±i aric? lim m{k)Y2{k) = -1. 

In particular, Yi and 12 have normalised poles at ±kc, and this allows us to derive the following 
linear and continuous inversion formula for A4. 

Lemma 9. Let Q be given with Q G L^(M) n BC^(M). Then there exists a unique Z G L^(M) such 
that 

M (z - l^'+*'> - ^'-'-'V. - + '^'-"'V.) = Q. (13) 

Moreover, Z depends linearly on Q and satisfies 

11^112 <c(||Q||2 + ||Q||l,oo) 

for some constant C independent of Q. 
Proof. The function Z with 

^(') Mk) 

is well-defined and continuously differentiable for k ^ ike- In view of Remark 8, I'Hospital's rule 
ensures that the limits limfc^_fc^ Z{k) and limfc_!>+fc^ Z{k) do exist, and combining this with the 
integrability properties of m and Q we find Z G L^(M). The inverse Fourier transform Z G L^(M) is 
therefore well-defined by Parseval's theorem, depends linearly on Q, and satisfies (13) by construction. 
With J := [—2kc, +2kc] we readily verify the estimates 

I|-^IIl2(R\J) < II"T-~"^IIl°°(]R\J)IIQIIl2(R\J) + (|Q(+fcc)| + |Q(-fcc)|) (||^i||l2(R\J) + II^2||l2(M\J)) 
<t?(l|Q||L2(M) + ||g||Loo(R)), 

and Taylor expanding both the numerator and the denominator of the right hand side in (14) at 
k = ±kc we get 

II-^IIl2(j) < C'll-^llco(j) < C'IIQIIci(j) • 

The desired estimate for \\Z\\2 now follows from ||Z||^2(ir) = II'^IIl2(m\j) ~^ W^Wl^^j)- Finally, Z is the 
unique solution in L^(M) since any other solution to (13) diflFers from Z by a linear combination of 
cos (kc ) and sin (kc-), see Remark 5. □ 
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Lemma 9 implies that the hnear operator M. admits a hnear and continuous inverse 
A^-^: jr-i(L2(M)nBC^(M)) ^ L2(M) span (Yi, Y2} , 



where J-~^ means inverse Fourier transform. The proof of Lemma 9 also reveals that can 
be extended to a larger space since one only needs that Q is continuously differentiable in some 
neighbourhood of zizkc- For our purpose, however, it is sufficient to assume that Q G BC^(M). We 
also mention that the constant C in Lemma 9, which is the Lipschitz constant of A4^^, is uniform in 
Co < c < ci but will grow with ci — >■ 1, due to the definition of Yi and Y2 and the properties of m. 



3.2 Solution operator to the cifRne subproblem 

We are now able to prove that the affine problem (10) 1 admits a solution operator 

£,: G eY ^ {S,r]) eXxR, 

where 

Y := |g G L°°(M) : suppG C [-1, 1]| . 

The existence of £ is a consequence of the following result. 

Lemma 10. For each G G Y there exists a unique (S, 77) G X x M such that 

MS = A'^G + r]. 
Moreover, S and 77 depend linearly on G and we have 



(15) 



(i) H<C^M'G||oo, 

(ii) \\S\\oo<Cm\\A^G\\oo, 
(Hi) ||S"||oo < Ca^II-^GIIoo , 

M ||S"||oo<C^||G||oo. 

for some constant Cm > independent of G. 

Proof. The function Q := J\?G satisfies suppQ C [—2, 2], and using 



Q{k) 



+ 



d 

dk 



Q{k) 



<C / {l + \x\)\Q{x)\ dx<C||Q||c 



as well as ||Q||2 = IIQII2, we easily verify that 

||Q||2 + ||Q||l,oo<C||Q||oo. 

By Lemma 9, the function S := M.~^A^G takes the form S = Z + /lYi + 72^2, where Z G 
and /i , /2 G M satisfy 

11^112 + l/ll + I/2I < CIIQIloo = C||^2^||oo . 

In particular, we have M.S = A?G and hence 

c^Z = J?Z - J?G + fxMYx + /2A^l2 • 



(16) 
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The functions MYi, MY2 are supported in [—1, +1], see Remark 8, and G € Y combined with 
Lemma 6 imphes that J^G vanishes outside of [—2, +2]. For |x| > 2 we therefore find 



c2 



.x+l/2 . s 2 ^ 



\Z{x)\ = \{A'z)iz)\ < {{AZ){s)) dsj 0, 



thanks to Holder's inequahty and since Lemma 6 imphes AZ G L^(M). By definition of M, Q, and 
S we also have 

c'^S = -A^G + A^S = -A^G + A^{Z + fiYi + /aYa) , (17) 
and Lemma 6 ensures that 

ll-^ '^||oo^||-^||2? Il^^lloo^ 11^ II 00 • 

Combining these estimates with (16) and (17), we arrive at 5 G L°°(M) with 

II'-' lloo < C||^^G|| 

Moreover, differentiating the first identity in (17) with respect to x, we get 

c^S' = V{-AG + AS) , c^S" = VV{-G + S) , 

where the discrete differential operator V is defined as VU = U(^- + ^) ~ " i)' Lemma 6. 
This implies 

11*5 lloo ^ C'll-^G'lloo 5 \\<S' lloo ^ CIICIloO 

thanks to ||^^G||oo < ||.AG||oo < ||G||oo and H^^Hoo < 11-511 00 ■ Since S does not belong to X, we 
now define 



\/2Tr \/2tt 
S{x) := S{x) - S{0) - h -\r\ ( cos {kcx) - l) - /2 ^7777^ sin {k^x) (18) 

m\kc) m'{kc) 

as well as 

and observe that S & X and (15) hold by construction. Moreover, S and rj depend linearly on G 
and the above estimates for /i, /i and 5 provide the desired estimates for both S and tj. Finally, 
the uniqueness of {S, rf) is a direct consequence of 5 G X and Lemma 5. □ 

Notice that the solution (5, rj) to (10)i is unique only in the space X x M and that further solution 
branches exists due to the nontrivial kernel functions oi M.. For instance, replacing (18) and (19) by 

S{x) := S - S{Q) , r/:=-(l-c2)5(0) 

we can define an operator 

Gg (5, ??) eXxM, X := (5 G W2'~(M) : 5(0) = O} , (20) 

which provides another solution to the affine problem (10)i. The corresponding corrector S", however, 
does in general not belong to X as it is oscillatory for both x — t- — 00 and x — t- +00. 

We emphasise that the three-parameter family of travelling waves R = Rq+S, which we construct 
below by fixed points arguments involving C, is - at least for sufficiently small S - independent of 
the details in the definition of C The reason is, roughly speaking, that changing C is equivalent to 
changing Rq, see the discussion at the end of Section 4. However, choosing X x M as image space for C 
provides more information on the resulting family of travelling waves: The existence of lima;^_|_oo S{x) 
reveals that for each c there exists one wave R = Rq + S that complies with the causality principle 
as it is non-oscillatory for x +00. 
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3.3 Properties of the nonlineeir operator Q 



In order to investigate the properties of the nonhnear superposition operator Q, we introduce a class 
of admissible perturbations S. More precisely, we say that S" G X is 6-admissible if there exist two 
number x_ < < a;+, which both depend on S and 6, such that 

(i) Ro{x±) + S{x±) = ±§ , 

(ii) Ro{x) + S{x) < -S for x < x- , 

(iii) Ro{x) + S{x) > +S for x > x^ , 



(iv) ^R'oiO) < R'oix) + S'{x) < 2R'q{0) for x- < x < x 



+ ) 



where Rq is the chosen wave for (5 = 0. Below we show that each sufficiently small ball in X consists 
entirely of (5-admissible functions, and this enables us to find travelling waves by the contraction 
mapping principle. 

We are now able to derive the second key argument for our fixed-point argument. 
Lemma 11. Let S eX be S-admissible and G = Q{S) as in (11). Then we have 

\\G\\oo<C, supp G Q[-CS, C6], [ Gix)dx < C{l + \\S"\\^)5^ 

for some constant C independent of S and S. 

Proof. The first assertion is a consequence of ||G||oo < 1 + C*- Since S is (5-admissible, we have 

supp G=[x-,x+], ±S = ± [Kix) + S'{x)) dx 

Jo 



with x± as above, and this implies 

1 



5 < \x±\ < 



i?'o(0) 



supp G C 



^o(O) 



-S,S] 



2i?'o(0) 
Using the Taylor estimate 

\R'oix) + S'ix) - R!o{0) - S'iO)\ < (||i?(;||oo + ||5"||oo) \x\ 



we also verify that 
x± T 



R'oiO) + S'iO) 
A direct computation now yields 

x+ 



^ ll-Rolloo + ||'S'"||oo ^ 11-^0 lloo "I" 11*^' 



'"I 



2 i?'o(0) + 5'(0) 



/ G{x) dx = [ ^ "^'siRoix) + S{x)) dx- f * sgn{Ro{x)) dx 

Jr Jx- Jx- 



^ dr 
-s 



z{r) 



(21) 



(22) 



(23) 



|a;+ + , (24) 



due to sgn(i?o(a;)) = sgn(x). Here, the function z with z[Ro{x) + S{x)) = Rq{x) + S'{x) for all 
X G x+j is well-defined since R + So is strictly increasing on [a:_, x+]. Thanks to (22), our 
assumption Is = J^g ^si"!") '^f = 0) ^.nd the estimate z{r), z{0) > ^Rq{0) we get 



^ dr 

-s 



z{r) 



-s 



*',(r) 



zlr 



z{0) 



dr 



< 



M \z{r)-z{0)\ ^ 
l*6(r)| ; s dr 



-5 



z{r)z{0) 



< CS 



{\x+\ + \X-\){\\R'^\\oo + \\S"\\oo) 
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and combining this with (21), (23) and (24) gives 
G{x) dx 



^\ , I , /-Ye/ I I I I I N ll-^olloo + ||'S'"||oo ^ ^e2 ll-^olloo + IliS"'! 

< |x_ + x+l + Go( + |a;+| j — . ^ < Go 



(25) 



By Proposition 2 (ii)c, i?o(0) bounded from below. Moreover, combining Proposition 2 (ii)a with 
the equation for i?Q, that is 

c^i?^; = Aii?o-Aisgn, 

we find a constant C, which depends only on cq and ci, such that ||i?o||oo !^ C . The second and third 
assertion are now direct consequences of these observations and the estimates (23) and (25). □ 



(26) 



Corollary 12. There exists a constant Cg, which is independent of S, such that 

MGIloo < CgS, M'Glloo <Cg{l + \\S"\\^)S'' , 

holds with G = G{S) for all 5-admissible S. 

Proof. Thanks to Lemma 11 and since A is the convolution with the characteristic function of the 
interval [— ^, +^], there exists a constant C such that 

\AG{x)\ < CS for \x±l\<CS, 

AG{x) = J^G{x)dx for \x\ < ^ - C6 , 
AG{x) = for \x\>l + C5, 

see Figure 5 for an illustration. The first bound in (26) is now a consequence of the trivial estimate 
|JjgG(a;)dx| < |suppG| ||G||oo < CS, whereas the second one follows from 



I [A'^G) (x) \<C6^+ f G{x) dx 
Jut 

and the refined estimate |/jgG(a;)dx| < C(l + ||5'||'4)'5^- 



for all X G 



□ 



jraph of G 



graph of AG 





























1.0 


1.0 1 









graph of A'^G 

be*') 




Figure 5: Properties of G = Q{S) for (5-admissible S. The shaded regions indicate intervals with 
length of order 0{5). 

In the general case Is ^ 0, one finds - due to /j^ G{x) dx = 2Is + 0{S^) - the weaker estimate 

||^^G||oo < C(l + ||5'^||)5. This bound is still sufficient to establish a fixed point argument, but 
provides only a corrector S of order 0(5). Recall, however, that Lemma 7 shows that shifting 
and changing Rq allows us to find correctors of order 0(5^) even in the case 1$ / 0. 

We finally derive continuity estimates for Q. 

Lemma 13. There exists a constant Cl independent of 5 such that 

\\A^G2 - A^GiWoo + \\AG2 - AGiWoo + <^||G2 - GilU < Cl^WS'^ - S[\\oo 
holds for all 6 -admissible correctors Si and S2 with G^ = Q{Se). 
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Proof. According to Lemma 11, there exists a constant C, such that Gi{x) = for all x with |a;| > C5. 
For |x| < C5, we use Taylor expansions for — 52 at a; = to find 

\S2{x) - Si{x)\ < \\S'-,-S[\\^\x\ , 

where we used that 52(0) — 51(0) = 0. Combining this estimate with the upper bounds for gives 

\G2{x) - Gi(x)| < y |52(x) - 5i(x)| < C\\S'^ - S[\\^ 
for all \x\ < CS, and this implies the desired estimate for \\G2 — Gi||oo- We also have 

11^2^2 - ^2^1 < \\AG2 - ^Gilloo < Isupp {G2 - Gi)\ \\G2 - Gi\\^ < G5\\G2 - Gllloo , 
and this completes the proof. □ 

3.4 Fixed point argument 

Now we have prepared all ingredients to prove that the operator 

r:=rsoCog 

admits a unique fixed point in the space 

:= G X : Halloo < (^ofJ^ , H^'Hoo < C'lcJ , IIS^Iloo < ^2! . 

Here, Ps denotes the projector on the first component, that means Vs{S, rj) = S, and the constants 
Gi are defined by 

C2 ■= Cm{^ + G>if), Ci := GmCq , Co := GMGg{l + G2) ■ 

Notice that any fixed point of T provides a solution to (10) and vice versa. 

Lemma 14. For all sufficiently small S, the operator T has a unique fixed point in Xs- 

Proof. Step 1: We first show that each 5 G X5 is (5-admissible provided that 5 is sufficiently small. 
According to Proposition 2, there exist positive constants tq, xq, and do such that 

|-Ro(3;)| > ''"0 for \x\ > , do < B'q{x) for |x| < xq , 

and combining the upper estimate for || -Roll 00 with the equation for Rq we find ||-Ro||oo ^ D2 for some 
constant D2. We now set 

" 2 "^^^ I 2f^' V ^' '7 ' ^= do^ 

and assume that 5 < 5q. For any x with |x| < a;^ < xo, we then estimate 

\R!q{x) + S'{x) - R!q{Q)\ < D2XS + G16 < (2^ + Giy< \do < iR'oiO) , 
and this gives iR'oiQ) < R'oix) + S'{x) < ^R'oiO). Moreover, xs < \x\ < xo implies 



\Ro{x) + S{x)\ > [ R'o{s)ds 
Jo 



2 

- US' lloo > {do — G16) \x\ > ^do ■ -^S = 6 , 

"0 



whereas for |x| > xq we find 

\Ro{x) + S{x)\ >ro- Gi6^ > ^ro > S . 
15 



Using 

X- := max{x : Ro{x) + S{x) < —6} , x+ := min{a; : Ro{x) + S{x) > +6} , 

we now easily verify that S is 5-admissible. 

Step 2: We next show that Ti^^s) C for all 5 < 6q. Since each 5 G X5 is (5-admissible, Corollary 
12 yields 

11^^(5) lloo < Cg5, ||^'^(5)||oo < Cg{l + C2)5^ 

and 11^(5') II 00 < 1 + C"^ holds by definition of Q and Assumption 1. Lemma 10 now provides 

||r(5)||oo < CMCg{l + C2)5'' = CoS\ 

l|r(5)'||oo < CMCgS = CiS, 

Iir(5)"||oo < Cm{i + c^) = C2, 

and hence T{S) G X^. 

Step 3: We equip X5 with the norm HS'H^ = II^Hoo + ||5"||oo + <^||'S"'||oo) which is, for fixed 5, 
equivalent to the standard norm. For given S'i,S'2 G X^, we now employ the estimates from Lemma 
10 and Lemma 13 for 5 = ^2 - Si and G = g{S2) - Q{Si). This gives 

||r(52) - r(-Si)||# < Cm {\\A^Q{S2) - AMs^)U + \\M{S2) - Ag{S^)U + ^11^(^2) - Q{S^)\\oo) 

< CmClS\\S2 — S'lWoo < CmCl^\\S2 — Si\\^ , 

and we conclude that T is contractive with respect to ||-||# provided that 6 < 1/{CmCl)- The claim 
is now a direct consequence of the Banach Fixed Point Theorem. □ 

The previous result implies the existence of a three-parameter family of waves that is parametrised 
by the speed c G [cq, ci] and by Rq, where Rq can be regarded as parameter in the two-dimensional 
L°°-kernel ot M. 

Proposition 15. Suppose that Is = for all 5. Then there exists Sq > with the following property: 
For any S < 60, each c G [cq, ci], and any Rq as in Proposition 2 there exists an 6-admissible 




SgXsH Il\R)® span |i, Yi - cos {kc-) , Y2 - 

such that R = Ro + S and solves the travelling wave equation (6) for some fi. In particular, we have 
-R(O) = 0, the limits 

lim (R(x) — a- cos (kcx) — P- sin (kcx)) and lim (R(x) — Ro(x)) 

x— >■— 00 s— >-+oo 

are well-defined for some constants a_, /3_ depending on c and Rq, and the estimates 

R{x)<-5 for x<-C5, R{x)>+6 for x>+CS (27) 
hold for som,e constant C > independent of c and Rq. 

Proof. For given c and Rq, Lemma 14 provides a unique fixed point S G of T, which solves 

MS = A^g{S) + r) 

for some r/ G M, and this implies that R = Rq + S is in fact a travelling wave. Moreover, by 
construction - see the proof of Lemma 10 - we also have 

\/2tt \/2Tr 

s = z + A + A (y, - ^ cos (*,)) + /. (y. - - (fc-)) 

for some constants /i, /2 and A and a function Z G L^(M) with Z{x) — t- as x — t- ±00. The claims on 
the asymptotic behaviour as x ^ ±00 now follow immediately since Rq has harmonic tails oscillations 
with wave number kc- Finally, the fixed point S is (5-admissible - see the proof of Lemma 14 - and 
this implies the validity of (27) due to < x+, —x- < C5. □ 



16 



Notice that Proposition 15 yields a genuine three-parameter family in the sense that different 
choices of the parameters c and Rq correspond to different tail oscillations for x +00 and hence 
to different waves R = Rq + S. This finishes the existence proof of Theorem 3. 

4 Uniqueness of phase transition waves 

In this section we establish the uniqueness result of the Theorem 3 by showing that the family 
provided by Proposition 17 contains all phase transition waves that have harmonic tails oscillations 
for X — > +00 and penetrate the spinodal region in a small interval only. 

Lemma 16. Let k > ^ be given and suppose that Is = for all S. Then there exists 6^ > such 
that the following statement holds for all < S < 5^: Let {Ri, /ii) and {R2, /X2) be two solutions to 
the travelling wave equation (6) with speed c G [cq, ci] such that 

Ri e w2'°°(M) , meR, Ri{0) = 

and 

Ri{x) < -S for x<-6'', Ri{x) > +5 for x > +6'' 
for both i = 1 and i = 2. Then, Ri and R2 are either identical or satisfy 

Ri{x) — R2{x) — a+( cos (fccx) — l) — /3+ sin {kcx) — 7+ y q 

for some constants 7+ and («+, 7^ (0, 0). 
Proof. For given Ri, R2, there exist constant /xi,/X2 € K such that 

M{R2 - Ri) = A'^G + fi2 - f^i , G := %{R2) - %{Ri) . 
By assumption and due to the bounds of we also find G{x) = for \x\ > 6'^ as well as 

\Gix)\<^\R2ix)-Ri{x)\<Cd''-^\\R!2-R[\\oo for \x\ < , 

and this implies 

\\AG\\^ < |suppG| ||G||oo < CS^^-^R'2 - R[\\oo ■ 
Moreover, Lemma 10 provides S eX as well as 77 e M such that 

MS = A'^G + ri, \\S'\\oo < C5^''-^\\R!2- R'lWoo- 

In particular, we have 

A^G = M{R2 - Ri - {I - c^y\l^2 - 1^1)) = m(^S - {1 - c^y\y 

Since the space of bounded kernel functions for M. is spanned by sin (kc-) and cos {kc-), we conclude 
that there exist constants a+ and /3+ such that 

R2{x) — Ri{x) = S{x) — a + a+(l — cos (kcx)) + /3+ sin {kcx) + 7+ , 

where a := lim2;^>+oo S{x) and 7_|_ := (l — c^) ^(//2 — IJ'i — v) + ~ (^+- case of a+ = /3+ = we 
therefore find 

1 1 r?' 7?' 1 1 1 1 Q' 1 1 ^ ^2k—1 \ \ jj/ p/ 1 1 

||-fi,2 -flj^lloo — 11'-' 1 1 00 _i L'O 1 1 -'''2 -'•'ill 00 ' 

and combining this with Ri (0) = R2 (0) we get R2 = Ri for all sufficiently small 5. □ 
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Proposition 17. Suppose that Is = for all 5 and that n with ^ < k <1 is fixed. Then there exists 
6k with < Sk < So such thai the following statement holds for all < 6 < 6^: Let R be a travelling 
waves with speed c G [cq, ci] such that the limit 

JunjRix) - Ro{x)) 

is well-defined for some Rq from Proposition 2 and such that 

R{x) < -S for x<-S'', R{x) > +6 for x > +6". 
Then R belongs to the family of waves provided by Proposition 15. 

Proof. Let i?o + 'S' be the travelling wave from Proposition 15. By construction, R — Rq — S converges 
as X ^ +00 and for all sufficiently small S we also have CS < S'^. Lemma 16 applied with Ri = R 
and R2 = Ro + S therefore implies R = Ro + S. □ 

With Proposition 15 and Proposition 17 we have established our existence and uniqueness result 
in the special case that Is = holds for for all S. The corresponding result for the general case is 
then provided by Lemma 7. 

We finally mention a particular consequence of our uniqueness result, namely that the family 
from Proposition 15 does not depend on the particular choice of the solution operator £. to the afHne 
problem (10)i. At a first glance, this might be surprising since the operator T and hence each fixed 
point surely depend on C We can, however, argue as follows. Suppose we would choose in the proof 
of Lemma 10 another reasonable solution operator C (for instance, the operator from (20) that does 
not involve any kernel function of M). Repeating all arguments from Section 3 we then find - for 
any given S, c, and Rq - a different corrector S G W^'°°(M). In general, this corrector S does not 
converge as a; ^ +00 but satisfies 

5(0) = 0, \\S\\oo<CS\ \\S'\\^<C6, ||5"||oo<C' 

for some constant C that is independent of c, Ro, and 5. Moreover, we also have 

S e L^(M) e span (l, Yi, Y2, cos (kc-), sin (kc-)} 

that means the tail oscillations of S for both x — > —00 and x — > +00 are again harmonic waves with 
wave number kc. Adding a suitable linear combination of 1 — cos (fee-) and sin(A;c-) to Ro we can 

construct another wave ^0 such that ^0 a-^d Rq + S have the same tails oscillations as x — > +00. 
This function ^0 is, at least for small 6, also a travelling wave for the unperturbed problem and 
hence among the family of waves provided by Proposition 2. We can therefore use ^0 instead of Ro 
in order to define the operator Q. Theorem 10, which relies on the oscillation-preserving operator C, 
then provides a corrector S that converges as x — t- +00, and from Lemma 16 we finally infer that 
Ro + S = Rq + S because both waves have, by construction, the same tail oscillations for x — t- +00. 
We therefore conclude, at least for small 6, that changing C does not alter the family of travelling 
waves but only its parametrisation by Rq. 

5 Kinetic relations 

We finally show that the kinetic relation does not change to order 0{5). To this end we denote by 
Rs a travelling wave solution to (2) as provided by Theorem 3. The corresponding configurational 
force, cf. [HSZ12], is then defined by := Te,j — Tf^^ with 

where the macroscopic strains fs,± on both sides of the interface can be computed from Rs via 

1 f^^ 

rs,±= hm - / Rs{±x)dx. 
L^oo L Jo 



18 



Lemma 18. Let Rs be a travelling wave from Theorem 3, and Rq the corresponding wave for 5 = 0. 
Then we have = Tq + 0{5^). 

Proof. By construction, we know that the only asymptotic contributions to the profile Rs are due to 
i?o — I 5 plus a small asymptotic corrector of order 0{5^) from span {l, Yi, 12}- This implies 

r5,± = ro,± - Is + 0{S'^). 

As fo,± and fs,± are both larger than S we know that 



Thus, we conclude 



and hence 



^'sirs,±) = rs,± T 1 = $o(ro,±) - Is + 0{S^) , 



Tf,5 = Tf,o - Is{ro,+ - ro ,-) + 0{6 ) . 

Moreover, we calculate 

Te,5 = ^'sir) dr = (r - ^'^(r) - *'o(r) + %{r)) dr 



$o(r) dr- (*5(r) - *o (0) dr 

•Jfx- 



rs,- -^rs,- 

= Mrs,+) - Mn,-) - 2Is = Urs,+ - 1)' - Un,- + 1)' - 2Is 
= i(fo,+ -Is- if - lifo,- -Is + if - 2Is + 0{S^) 
= i(fo,+ - 1)' - i(fo,- + 1)' - Is{fo,+ - 1 - fo,- - 1) - 2Is + 0{S^) 
= Te,o - Is{fo,+ - fo,-) + 0{6^) . 

Subtracting both results gives T5 = Tq + 0{6^), the desired result. □ 
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